Assessing population genetic structure and variability with RAPD data: Application to Vaccinium macrocarpon (American Cranberry)

A method for estimating and comparing population genetic variation using random amplified polymorphic DNA (RAPD) profiling is presented. An analysis of molecular variance (AMOVA) is extended to accomodate phenotypic molecular data in diploid populations in Hardy-Weinberg equilibrium or with an assumed degree of selfing. We present a two step strategy: 1) Estimate RAPD site frequencies without preliminary assumptions on the unknown population structure, then perform significance testing for population substructuring. 2) If population structure is evident from the first step, use this data to calculate better estimates for RAPD site frequencies and sub-population variance components. A nonparametric test for the homogeneity of molecular variance (HOMOVA) is also presented. This test was designed to statistically test for differences in intrapopulational molecular variances (heteroscedasticity among populations). These theoretical developments are applied to a RAPD data set in Vaccinium macrocarpon (American cranberry) using small sample sizes, where a gradient of molecular diversity is found between central and marginal populations. The AMOVA and HOMOVA methods provide flexible population analysis tools when using data from RAPD or other DNA methods that provide many polymorphic markers with or without direct allelic data.     

Potential pitfalls of the PTV concept in dose-to-medium planning optimization

In typical treatment planning of 3D IMRT, the incident energy fluence is optimized to achieve a homogeneous dose distribution to the PTV. The PTV includes the tumour but also healthy tissues that may have a different dose response for the same incident energy fluence, like bony structures included in the PTV (mandibles in head and neck tumours or femoral bones in sarcomas). Dose to medium optimization compensates for this heterogeneous response, leading to a non-homogeneous energy fluence in the PTV and a non-homogeneous dose in the CTV in the presence of geometric errors. We illustrate qualitatively this statement in a cylindrical geometry where the PTV includes a CTV (7 cm diameter) made of water surrounded by ICRU compact bone (1.2 cm thickness); such configuration was chosen to exaggerate the aforementioned effect. Optimization was performed assuming dose equals photon energy fluence times mass energy absorption coefficient. Bone has a 4% lower dose response in a 6 MV flattening filter free spectrum. After optimization either in medium or assuming everything as water composition, the geometry was shifted by 1.2 cm and dose recomputed. As expected, compensating for the under-response of the bone material during optimization in medium leads to an overdosage of the CTV when patient geometric errors are taken into account. Optimization in dose assuming everything as water composition leads to a uniform coverage. Robust optimization or forcing a uniform atomic composition in the PTV margin may resolve this incompatibility between the PTV concept and dose to medium optimization.     

On the choice of the optimal periodic operation for a continuous fermentation process

In this contribution we investigate the impact of the forcing waveform on the productivity of a continuous bioreactor governed by an unstructured, nonlinear kinetic model. The (periodic) forcing is applied on the substrate concentration in the feed. To this end, some alternative waveforms commonly encountered in practice are evaluated and their performance is compared. An analytical/numerical approach is used. The preliminary analytical step is based on the π‐criterion that gives useful information for small amplitudes. The extension to larger amplitudes, when significant improvements are expected, is then performed through a continuation‐optimization procedure. It is found that the choice of the specific waveform has an impact on the performance of the process and there is no unique best forcing for any process condition, but its choice depends on the operating parameters and the forcing amplitude and frequency values. Further, the influence of the waveform functions on the wash‐out conditions are extensively examined. The analysis shows that all the waveforms examined in this work may lead to significant enlargement of the nontrivial regime with respect to a steady state operation. In particular, square‐wave forcing leads in practice to the extinction of the wash‐out conditions for any feed substrate concentration and for a well defined choice of the forcing parameters. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Construction of Selection Indexes with Restrictions

In this paper the problem of unistage selection with inequality constraints is formulated. If the predictor and criterion variables are all normally distributed, this problem can be written as a convex programming problem, with a linear objective function and with linear constraints and a quadratic constraint. Using the duality theory, for convex nonlinear programming it is proved, that its dual problem can be transformed into a convex minimization problem with non-negativity conditions. Good computational methods are known for solving this problem. By the help of the dual problem sufficient conditions for a solution of the original primal problem are derived and illustrated by an example of practical interest.     

Hydraulic acclimation in a Mediterranean oak subjected to permanent throughfall exclusion results in increased stem hydraulic capacitance

Stem water storage capacity and hydraulic capacitance (C_S) play a crucial role in tree survival under drought-stress. To investigate whether C_S adjusts to increasing water deficit, variation in stem water content (StWC) was monitored in vivo for 2 years and related to periodical measurements of tree water potential in Mediterranean Quercus ilex trees subjected either to permanent throughfall exclusion (TE) or to control conditions. Seasonal reductions in StWC were larger in TE trees relative to control ones, resulting in greater seasonal C_S (154 and 80 kg m⁻³ MPa⁻¹, respectively), but only during the first phase of the desorption curve, when predawn water potential was above −1.1 MPa. Below this point, C_S decreased substantially and did not differ between treatments (<20 kg m⁻³ MPa⁻¹). The allometric relationship between tree diameter and sapwood area, measured via electrical resistivity tomography, was not affected by TE. Our results suggest that (a) C_S response to water deficit in the drought-tolerant Q. ilex might be more important to optimize carbon gain during well-hydrated periods than to prevent drought-induced embolism formation during severe drought stress, and (b) enhanced C_S during early summer does not result from proportional increases in sapwood volume, but mostly from increased elastic water.     

On the parametrization of lateral dose profiles in proton radiation therapy

PurposeThe accurate evaluation of the lateral dose profile is an important issue in the field of proton radiation therapy. The beam spread, due to Multiple Coulomb Scattering (MCS), is described by the Molière's theory. To take into account also the contribution of nuclear interactions, modern Treatment Planning Systems (TPSs) generally approximate the dose profiles by a sum of Gaussian functions. In this paper we have compared different parametrizations for the lateral dose profile of protons in water for therapeutical energies: the goal is to improve the performances of the actual treatment planning.MethodsWe have simulated typical dose profiles at the CNAO (Centro Nazionale di Adroterapia Oncologica) beamline with the FLUKA code and validated them with data taken at CNAO considering different energies and depths. We then performed best fits of the lateral dose profiles for different functions using ROOT and MINUIT.ResultsThe accuracy of the best fits was analyzed by evaluating the reduced χ², the number of free parameters of the functions and the calculation time. The best results were obtained with the triple Gaussian and double Gaussian Lorentz–Cauchy functions which have 6 parameters, but good results were also obtained with the so called Gauss–Rutherford function which has only 4 parameters.ConclusionsThe comparison of the studied functions with accurate and validated Monte Carlo calculations and with experimental data from CNAO lead us to propose an original parametrization, the Gauss–Rutherford function, to describe the lateral dose profiles of proton beams.     

A genetic linkage map of Theobroma cacao L.

A linkage map of the cocoa genome comprising 193 loci has been constructed. These loci consist of 5 isozymes, 101 cDNA/RFLPs, 4 loci from genes of known function, 55 genomic DNA/RFLPs and 28 RAPDs. A population of 100 individuals derived from a cross between two heterozygous genotypes was used. Segregation analyses were performed with the JoinMap program. Ten linkage groups, which putatively correspond to the ten gametic chromosomes of cocoa, were identified. The map covers a total length of 759 cM with a 3.9 cM average distance between 2 markers. A small fraction (9%) of the markers deviated significantly from the expected Mendelian ratios.     

基于多尺度协同的长沙市生态网络构建与层级优化

城市化背景下人类与自然环境的矛盾呈现出多尺度、层级化特征，而传统生态网络的构建方式较少考虑不同尺度下生态要素的关系，无法从区域落实到中心城区，难以形成系统性的解决方案。研究在综合梳理各尺度生态网络构建方法的基础上，以长沙市为例，基于形态学空间格局分析（Morphological Spatial Pattern Analysis，MSPA）、景观连通性原理和生态斑块重要性评价识别生态源地，并通过多层级生态阻力面的确定，综合运用最小费用路径（Least-cost path method，LCP）、电路理论、层级传导理论、尺度嵌套等方法对市域、都市区、中心城区的生态网络进行了协同构建和层级优化，最后基于不同尺度生态网络的特点应用并落实到多层级的国土空间规划体系中。研究结果表明：（1）长沙市域生态网络和都市区生态网络具有较好的层级嵌套特征；共识别两尺度生态叠合源地14个、生态叠合廊道15条，主要通过中心城区内的湘江、浏阳河和捞刀河部分河段与外围生态绿圈相衔接，形成"外环内楔"的空间格局。（2）确定市域重要廊道、市域潜在廊道、生态叠合廊道、都市区重要廊道、都市区潜在廊道的核心保护面积共501.14 km²，并提取位于生态廊道核心保护区范围内的生态夹点和生态障碍点，以进一步落实生态保护修复策略。（3）得到具有重要生态连通功能的中心城区生态绿道长度441.2 km，生态修复单元56个，并结合生态阻力值划分为5级进行针对性修复。（4）基于不同尺度生态网络的衔接、嵌套，最终构建"市域总体生态安全格局-都市区城市生态空间发展格局-以城市绿道为基础的中心城区生态修复单元"，并与不同层级的国土空间规划体系相对应。研究结果将为以大城市为中心的跨尺度生态系统修复和生态安全格局构建提供科学参考。